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Context. The X-ray burster 4U 1850-087, located in the Galactic globular cluster NGC 6712, is an ultracompact binary (orbital period ~21 
minutes), likely harbouring a degenerate companion. 

Aims. The source has been observed at soft y-rays several times with the INTEGRAL satellite, during the monitoring of the Galactic plane, with 
an unprecedented exposure time. We analysed all available INTEGRAL observations, with the main aim of studying the long-term behaviour of 
this Galactic bulge X-ray burster. 

Methods. The spectral results are based on the systematic analysis of all INTEGRAL observations covering the source position performed 
between March 2003 and November 2005. 

Results. The source X-ray emission is hard and is observed, for the first time, up to 100 keV. A broad-band spectrum obtained combining 
the INTEGRAL spectrum together with a quasi-simultaneous XMM-Newton observation performed in September 2003 is well modeled with 
a disk-blackbody emission (with an inner disk temperature of ~0.8 keV) together with a power-law (with a photon index, F, of 2.1). The 
2-100 keV luminosity is 1.5x10^^ erg s"' (assuming a distance of 6.8 kpc). 

Conclusions. INTEGRAL observations reveal for the first time that this X-ray burster displays a very hard X-ray spectrum, with a cut-off at 
energies higher than 100 keV, and that the source spends most of the time in this low luminosity and hard state. Indeed, a previous BeppoSAX 
observation in April 1997 observed high energy emission from 4U 1850-087 only up to 50 keV. 
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1. Introduction 

4U 1850-087 (Swank et al. 1976) is an X-ray burster located 
in the galactic globular cluster NGC 6712. The likely optical 
counterpart displays an UV modulation (Anderson et al. 19931) 
with a period of 20.6 minutes (interpreted as the orbital pe- 
riod), implying a degenerate companion of O.O4M0 (Homer et 



al. [T996l l. 

The first broad-band spectrum from this source was ob- 
tained with BeppoSAX in the energy range 0.3-50 keV during a 
pointing performed in 1997 (Sidoli et al. '2001V The spectrum 
was well described with a disk-blackbody and Comptonized 
continuum with A^h = 3.9 x 10^' cm"^, an inner disk temper- 
ature, kTin, of ~0.6 keV, an inner projected radius of ~5 km 
(for an assumed NGC 6712 distance of 6.8 kpc, HaiTis 1996 1, 
a temperature, kTo, of the input "seed" photons, of 0.8;^^ j ^eV, 
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an electron temperature, kT^, of 70 keV, and an optical 
depth, T, of 1.7. The estimated 0.1-100 keV luminosity was 
1.9x10^^ erg s-i (at 6.8 kpc). 

Other X-ray observations are reported for 4U 1850-087 
limited to the energy range below 10 keV (ASCA, Juett et 
al.IMm XMM-Newton, Sidoli et al. 110051 Chandra, Juett & 
Chakrabartv l2005t . 

We report here the first detection of 4U 1850-087 above 
50 keV with the INTEGRAL satellite, obtained with an un- 
precedented exposure time collected during the monitoring of 
the Galactic plane. 

2. Observations and Results 

The ESA INTEGRAL gamma-ray observatory, launched in 
October 2002, carries three co-aligned coded mask telescopes: 
the imager IBIS (Ubertini et al. 2003), which allows high an- 
gular resolution imaging over a large field of view (29° x 29°) 
in the energy range 15keV-10MeV, the spectrometer SPI 
(Vedrenne et al. 2003; 20 keV-8MeV) and the X-ray moni- 
tor JEM-X (Lund et al. 2003; 3-35 keV). IBIS is composed of 
a low-energy CdTe detector (ISGRI; Lebrun et al. 2003), sensi- 
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tive in the energy range from 15 keV to 1 MeV, and a Csl detec- 
tor (PICsIT; Labanti et al. 2003), designed for optimal perfor- 
mance at 5 11 keV, and sensitive in the 175 keV-10 MeV energy 
range. 

We analyzed all public and Core Program IBIS observa- 
tions pointed within 10° of the source. This resulted in 909 indi- 
vidual pointings (Science Windows, SWs) performed between 
March 2003 and November 2005. All the data have been pro- 
cessed using version 5.1 of the OSA INTEGRAL analysis soft- 
ware, and analysed with the corresponding rebinned response 
matrices. 

A mosaic image of all 909 pointings revealed the sources 
that were active in this field during the period of interest. The 
list of all the detected sources (detection significance > 5cr) was 
then used in the spectral step, where spectra for all the detected 
sources in the field of view were simultaneously extracted with 
the standard spectral extraction method. 
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Fig. 1. RXTE/ASM Ughtcurve (day average) of 4U 1850-087. 
The small solid circles in the upper region show the times of 
the INTEGRAL observations analysed here. 



Table 1. Summary of all INTEGRAL observations of 
4U 1850-087 analysed here. The observations have been 
grouped together in 4 data-sets for brevity. The 3rd and 4th 
columns list the Start and Stop Time of the four groups of ob- 
servations, the fifth column reports the number of SWs in each 
data-set. Data-set n. 4 includes Core Program data still not pub- 
licly available. 



(with a null hypothesis probability in the range 0.1-2.4%) 
which can be explained by the spectral distortions present when 
extracting IBIS/ISGRI spectra spanning very different mis- 
sion times (Lubinski et al. 2005), especially for faint sources, 
and by possible intrinsic source spectral changes on smaller 
timescales, which cannot be a priori excluded. Nevertheless, 
the global spectral shape is rather constant on timescales of 
months-years. 



Data 


Temporal 


Start Time 


End Time 


Num of 


set 


window 


(MJD) 


(MJD) 


SWs 


1 


Mar 2003-May 2003 


52708.7 


52772.0 


258 


2 


Sep 2003-Nov 2003 


52910.9 


52963.4 


210 


3 


Mar 2004-May 2004 


53075.1 


53126.5 


84 


4 


Aug 2004-Nov 2005 


53238.1 


53684.0 


357 



Table 2. Results of the spectral analysis of the four temporal- 
selected IBIS/ISGRI spectra. Data-set numbers are the same as 
in Tabled The 2nd column reports the net exposure time of 
each spectrum. The 3rd column lists the IBIS/ISGRI average 
rates in the energy range 20-100 keV. Fluxes are in units of 
10-"' erg cm-2 s-i(20-100 keV). 



4U 1850-087 is a relatively faint source, therefore a mean- 
ingful spectral analysis needs to be performed adding together 
several observations. For this reason, we grouped the pointings 
in four data-sets (listed in Table Q from which we extracted 
four different IBIS/ISGRI spectra. 

The 2-10 keV flux of 4U 1850-087 shows a brightening 
event peaking on 25 February 2004, and lasting ~10 days, 
as visible from the ASM/RXTE lightcurve reported in Fig.^ 
Only part of its decay is covered by the INTEGRAL observa- 
tions of data-set 3. 

When we are averaging the spectra of 4U 1850-087 using 
data that span over different mission times, we use a system- 
atics of 5% in the spectral fit process, to take into account 
for modifications in the instrument evolution (Lubinski et al., 
2005). 

The best-fit spectral results, using a single power-law are 
reported in Tab.|2]for the four data-sets. The fits display a high 



Data- set 


Exp. 


IBIS/ISGRI rate 


Photon Index 


Flux 


;r'/dof 




(ks) 


(s-') 








1 


372 


1.25±0.05 


2.5+0.2 


0.8 


21.6/8 


2 


383 


2.36±0.05 


2.1±0.1 


1.8 


18.8/7 


3 


114 


1.42+0.09 


2.2±0.3 


1.1 


11.6/7 


4 


490 


1.59±0.05 


2.2+0.1 


1.2 


16.1/7 



We tried also to extract JEM-X spectra, but the smaller field 
of view, the faintness of the source, together with some instru- 
mental issues (e.g. JEM-X switched-off in several occasions), 
severely reduced the number of useful observations. Thus data 
could be used only for the integrated spectrum covering the 
data-set 3 and part of the data-set 4. The JEM-X spectra corre- 
spond to 76.7 ks net exposure time for JEM-Xl and 55.5 ks for 
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Fig. 2. 4U 1850-087 combined JEM-X and IBISASGRI simul- 
taneous spectrum, fit with a single power-law (photon index of 
2.2+0.2). In the lower panel, the residuals in units of standard 
deviations are shown. 

JEM-X2, with a source count rate of 0.23+0.02 and 0.22+0.03, 
respectively in the two JEM-X units (5-20 keV). 

The combined JEM-X plus IBIS/ISGRI simultaneous spec- 
trum (taken from the data-set 3) is shown in Fig. 12 fitted with 
an absorbed power-law, with absorbing column density fixed at 
2.5x 10^' cm"^ (the expected interstellar absorption towards the 
globular cluster NGC 6712). The fit is a good deconvolution of 
the spectrum (x^ ldof-21 SjlJ) and results in a photon index of 
2.2+0.2 and in a 5-100 keV flux corrected for the absorption 
of (2.4+0.3)xl0"'° erg cm^^ s"' (based on the IBIS/ISGRI re- 
sponse matrix). Other simple models, like bremsstrahlung or a 
cut-off power-law, never result in a better fit. 

4U 1850-087 has been observed with XMM-Newton on 27 
September 2003 (see Sidoli et al. 2005), i.e. within the period 
covered by INTEGRAL data-set 2, with the main aim to study 
the low energy absorption intrinsic to the source. Indeed the 
XMM-Newton spectrum below ~2 keV is complex and we re- 
fer to Sidoli et al. 2005 for its detailed description (and for 
the XMM-Newton data reduction). EPIC PN operated in Small 
Window mode, and with a net exposure time of 8. 1 ks. A com- 
bined spectral analysis EPIC/IBIS interestingly allows us to ex- 
tend the spectral study of this source in the soft X-rays, be- 
low 5 keV, and to significantly refine the spectral parameters. 
During the XMM-Newton observation, only EPIC PN data did 
not suffer from pile-up problems, thus we use this spectrum 
(1.7-12 keV, see Sidoli et al. 2005 for the details) to perform 
a broad-band analysis of the quasi-simultaneous EPIC-IBIS 
observations. 

Among the different models considered by Sidoli et al. 
(2005) when fitting the XMM-Newton spectrum, we tried their 
best-fit, a disk-blackbody emission together with a power-law. 
This gives already a very nice result (;t'^/dof= 147.4/203), with 
the spectral parameters reported in Table |3l and the spectrum 
shown in Fig.|3l The normalization constant for the IBIS/ISGRI 
spectrum relative to PN (constant factor fixed at 1), was 1 
We note that this normalization factor between the two instru- 
ments is in agreement with that reported for the Crab spectrum 
(Kirsch et al. l2005l . The broad band flux (2-100 keV, corrected 



for the absorption) is 2.8x10 erg cm ^ s ' (based on the 
EPIC PN response matrix). Fixing the column density to the 
best-fit (5.7x10^' cm"^) found by Sidoli et al. (2005), does not 
change the spectral results. 

In order to obtain physical information, we then adopted 
the best-fit model used for the BeppoSAX 4U 1850-087 spec- 
trum (Sidoli et al. 2001), a disk-blackbody plus a Comptonized 
emission (compTT model in xspec; Titarchuk ll994> . Since the 
data extend only down to 1.7 keV, we linked the tempera- 
ture of the seed photons, kTo, to the inner disk temperature, 
kTin. This is also justified by the empirical finding that in 
ultra-compact binaries the two temperatures are similar, within 
the uncertainties (Sidoli et al. 2001). We obtain a good fit 
(;t'-/dof=156. 3/203), with the parameters reported in Table |3 
To obtain better constraints to the spectral paramenters, in 
this case the absorbing column density has been fixed to the 
NGC 6712 interstellar absorption value. With this model, the 
normalization constant of the IBIS/ISGRI spectrum with re- 
spect to EPIC PN lies in the range 1.5-1.9, again in agreement 
with the standard value (Kirsch et al. 2005). Adopting a cut-off 
power-law instead of a power-law in the two-component model 
OifVdof= 172.2/202), the cut-off energy resulted in Ec> 1 10 keV 
(90% confidence level), with a best-fit photon index of 1 .9+0. 1 . 



3. Discussion and Conclusions 

Previously, the source broad-band spectrum was observed only 
with BeppoSAX in April 1997. At that time, the PDS non- 
imaging instrument, covering the range from 20 to 200 keV, 
detected the source only up to 50 keV (Sidoli et al. 2001). 

We report here for the first time the discovery of hard (50- 
100 keV) X-ray emission from the X-ray burster 4U 1 850-087 
and a long-term study of its X-ray spectral behaviour. It is now 
possible to directly measure the high energy luminosities of 
this source: they are ~1.3xlO-'^ erg s"' and 7x10-'^ erg s"' re- 
spectively in the 1-20 keV and 20-200 keV energy ranges (at 
6.8 kpc). This impHes that 4U 1850-087 falls into the so-called 
"burster box" (BaiTet et al., 2000). 

In order to allow a proper comparison with the INTEGRAL 
spectrum we re-analysed the BeppoSAX observation (April 
1997; refer to Sidoli et al. 2001 for the details of the data re- 
duction) adopting the same model used here, a disk-blackbody 
emission together with a power-law. This resulted in the fol- 
lowing parameters (x^/dof-llQA/lS^): an absorbing column 
density of (0.46+0.03)xl0^^ cm"^, an inner disk blackbody 
temperature, kTin, of 0.66+0.03 keV, and a powerlaw photon 
index of 1.96+0.06. Thus, there is a good agreement with the 
INTEGRAL spectroscopy. 

The new observations with INTEGRAL allow us to put 
much more stringent limits on the presence of a high energy 
cut-off. In Fig. 13 we compare the confidence contour levels 
for the high energy cut-off and the power-law photon index, 
obtained with the EPIC-IBIS joint spectrum (solid contours) 
and with BeppoSAX (dashed contours). During the BeppoSAX 
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Fig. 3. Broad-band 4U 1850-087 EPIC/PN spectrum (obtained in September 2003) combined with the quasi-simultaneous 
IBIS/ISGRI (data-set 2). On the left, the count spectra are displayed, together with residuals in units of standard deviations, 
when fitted with a disk-blackbody plus a power-law. On the right, the corresponding 1.7-100 keV energy spectrum is shown. 

Table 3. Results of the spectral analysis of the broad-band quasi-simultaneous spectrum (EPIC/PN-IBIS/ISGRI). Fluxes (2-100 
keV) are corrected for the absorption, are in units of 10"^ erg cm"^ s and are based on the EPIC pn response matrix. The 
assumed distance is 6.8 kpc (Harris et al. 1996). The compTT model assumes a spherical geometry. / is the inclination angle of 
the disk. 



Model 


Nh 

(10^2 cm-2) 


Parameters 


Flux 


xVdof 


DISKBB+POW 


0.4±0.2 


kTi„=0.8±0.1 keV 
ri„x(co,5(0)'" = 1.7!iikm 


2.8±0.1 


147.4/203 


diskbb+compTT 


0.25 fixed 


kTe>40 keV 
kTo=0.67±0.05 keV 
T=2.2+0.5 
kTi„=kTo 

ri„x(coi(!))°^=4±l km 


2.9±0.1 


156.3/203 



observation a lower limit to the high energy cut-off could be 
placed at Ec>60 keV (90% level), while now with INTEGRAL 
we can interestingly shift it towards much higher energies, 
with Ec>100 keV. This allows us to include 4U 1850-087 
among the few tens of low-luminosity low-mass X-ray bina- 
ries (LMXRBs) for which the spectrum has been observed up 
to 100 keV (Barret et al. 2000, Di Salvo & Stella 2002, Bazzano 
et al. 2006). More interestingly, 4U 1850-087 is now among 
the hardest type 1 X-ray bursters in our Galaxy. 

The INTEGRAL spectrum displays properties which fits 
well into the classification of low luminosity (~0.01 LEdd) 
weakly magnetized neutron stars with "hard spectra" (Barret 
et al. 2001): a broad-band spectrum extending up to 100 keV, 
well described by a soft disk-blackbody emission together with 
a hard Comptonized component, with a optical depth of the 
Comptonizing corona of ~2, seed photons temperature be- 
low 1 keV, a high electron temperature, above 40 keV. A 
few LMXRBs appear to spend most of the time in this "hard 
state" (Di Salvo & Stella 2002), like 4U 0614+091 (Piraino 
et al. 1999), which also displays spectral parameters similar 
to 4U 1850-087, with a photon index in the range 2.3-2.4. 



Also for 4U 1850-087 the four diff^erent INTEGRAL spectra 
suggest that the source spend most of its lifetime in this spec- 
tral state (on timescales of months or years, if compared with 
BeppoSAX). 

The very hard emission from 4U 1850-087 is clearly evi- 
dent also when compared (see Fig.|5} to a typical spectrum of 
a Z-source, GX 349+2, reported in an INTEGRAL study of the 
spectral behaviour of the persistent LMXRBs hosting a neutron 
star (Paizis et al. 2006). 
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